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We study spin chirality for a spin-^ zigzag XY chain at low temperature by applying a low- 
temperature density-matrix renormalization group technique. We calculate temperature dependence 
of dynamical and static correlations of the chirality. In a chiral phase, chiral long-range order at zero 
temperature disappears at finite temperature, consistent with the fact that there is no long-range 
order at finite temperature in one-dimensional systems with short-range interactions. In a dimer 
phase next to the chiral phase, we find an enhancement of static chiral correlation as well as spin 
correlation with increasing temperature. The enhancement corresponds to the increase in spectral 
weight inside a gap in the dynamical chiral-correlation function. This temperature-induced chiral 
correlation is a demonstration of the presence of a chiral state in excited states. 
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I. INTRODUCTION 

Recently low-dimensional quantum spin systems have 
been attracting much attention, since there are many in- 
teresting phenomena associated with magnetic frustra- 
tion. A spin-i zigzag XX Z chain is a typical one- 
dimensional (ID) quantum spin system with frustra- 
tion caused by competition between the first and second 
nearest-neighbor exchange interactions. Recent studies 
have shown that a vector chiral long-range order (LRO) 
(Ref. 1) occurs under certain conditionsi^"— More detailed 
studies of the model have been done very recently by 
searching for the vector chiral order in the whole param- 
eter region!^ It has also been shown that a magnetic field 
induces the vector chiral order even for a spin-i zigzag 
Heisenberg chain<^"— Therefore, the vector chirality can 
be one of important concepts characterizing ID spin-i 
frustrated magnets. The vector spin chirality is believed 
to be crucial for understanding the physics of ID mul- 
tiferroic cuprate o^'^° such as LiCu202 (Refs. 11-13) and 
LiCuV04 (Refs. 14-16). 

The ground-state properties of the spin-^ zigzag XX Z 
chain have extensively examined numerically^ and an- 
alytically.— However, physical properties at finite tem- 
peratures have not been studied as far as we know. Even 
for a simpler case without anisotropic z component, i.e., 
zigzag XY chain, there is no work on temperature de- 
pendence of physical quantities. In particular, it is inter- 
esting to clarify how the chiral properties change at finite 
temperature. 

In the chiral phase of the spin-i zigzag XXZ chain, 
it is known that time-reversal symmetry and global U{1) 
symmetry are preserved but parity symmetry is broken. 
It is interesting to clarify how finite temperature affects 
on such a chiral phase and chiral properties in the model. 

In this paper, we investigate finite-temperature prop- 
erties of the spin-i zigzag XY chain that is the simplest 
frustrated quantum spin model. The model shows vec- 
tor chiral LRO at zero temperature in some parameter 



regionsii It is known that the chiral ordered phase is lo- 
cated next to a dimer phasei^ In order to clarify the ef- 
fect of temperature on the chirality in both the chiral 
and dimer phases, we employ a recently developed low- 
temperature density matrix renormalization group (LT- 
DMRG) techniquei^ and calculate temperature depen- 
dence of dynamical and static correlations of the chirality. 
In the chiral phase, we find that the chiral LRO disap- 
pears at finite temperature as expected form the fact that 
there is no LRO at finite temperature in ID systems with 
short-range interactions. In the dimer phase, static chiral 
correlation as well as spin correlation is enhanced with 
increasing temperature, followed by subsequent suppres- 
sion with further increase in temperature. The maximum 
of the correlations appears at a temperature where spec- 
tral weight inside a gap in the dynamical chiral correla- 
tion function increases. This temperature-induced chiral 
correlation is a demonstration of the presence of a chiral 
state in excited states. 

This paper is organized as follows. The spin-i zigzag 
XY model is introduced in Sec. II. In Sec. Ill, LT- 
DMRG is briefly explained and quantities examined 
arc introduced. The dynamical chiral correlation func- 
tions at both zero and finite temperatures are shown in 
Sec. IV. A. The static correlation functions are presented 
in Sec. IV. B. Conclusions are summarized in Sec. V. 



II. MODEL 

The Hamiltonian of the spin-^ zigzag XY chain is 
given by 
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^ ^'^•^p'^ {^^l^l+p + Sf^l+p) y (1) 
p=l / 

where Sf^^^ is the x{y) component of spin operator at site 
I, and Jp is the nearest {p = 1) and next nearest (p = 2) 
neighboring spin coupling constant. In the present paper. 
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we take antiferromagnetic spin couplings, i.e., Ji > and 
J2 > 0, which lead to a frustrated spin system. At zero 
temperature, three phases have been identifieds^ gaplcss 
chiral phase (aj = J1/J2 ^ 0.81), gapped dimer phase 
(0.81 ^ ctj 3.1), and Tomonaga-Luttinger liquid phase 
{aj < 3.1). In the chiral phase, vector chirality 



Ki — S[Sf_^_^ — SfSf^j^ (2) 
shows a finite expectation value that is evidence of LRO. 



III. METHOD 

In the present study, we examine effects of both quan- 
tum and thermal fluctuations on the frustrated spin-i 
zigzag XY chain. However, in such a frustrated system, 
it is difficult to employ standard Jordan- Wigner trans- 
formation, because phase factor corresponding to a next- 
nearest neighbor interaction remains after the transfor- 
mation. Thus, we use the DMRG methodic to analyze 
properties of this frustrated system. 

In the standard DMRG method at zero temperature, 
we consider the ground state as a target state. At finite 
temperatures, on the other hand, we employ a state |^) as 
a target state. The state |^) is defined as a product of the 
Boltzmann factor and a normalized state |^) randomly 
generated)^ namely. 



(3) 



where /3 = 1/T with temperature T. (Hereafter we take 
the Boltzmann constant fee = 1 and h = 1.) Since an 
arbitrary state can be expanded with a complete set of 
basis the target state |^) is rewritten by 



TV 
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(4) 



where N is the total number of the basis and = (CilO 
is a normalized random number. In our procedure, is 
constructed by the basis of a total system in DMRG pro- 
cedure called superblock. |^) is obtained by the regulated 
polynomial expansion of the Boltzmann factorJ^ 

We next explain how to calculate thermodynamical 
quantities. \Ci) can also be expanded by eigenstate |e„) 
with eigenenergy e„. |<^) thus reads 



N N 



(5) 



2—1 71—1 



with bi^n = (eniC 



We consider a normalized expecta- 
tion value of an arbitrary operator A in terms of |^), 



(6) 



By using Eq. ([5]), the denominator is given by 



(7) 



If the coefficient is randomly generated from a rect- 
angular distribution and a plenty of the random num- 
ber set is averaged, a square of converges on a non- 
zero constant {rf)r = C, where (• • • denotes random 
sampling and averaging. A product {rirj)r converges on 
zero because positive and negative values appear with 
same probability during averaging. After random sam- 
pling and averaging, the denominator is expected to be a 
product of the constant C and the partition function Z, 



(8) 



because of the normalized condition J^i^in ^ 1- The 
numerator of ((A)) in Eq. ([6]) reads 



(9) 



where 



and 



Sn,' 



-/3(£„+e™)/2/ 



n I ^ I CjTj ) 



(10) 



(11) 



Similar to the denominator, after random sampling and 
averaging, the bracket [• • • ] in Eq. (fTO|) converges on the 
constant C. On the other hand, the bracket in Eq. pT|) 
converges on zero because (rf)r = C, {rirj)r — 0, and 
X^i bi,nbi^m = Sn,m- Thc numerator of {{A)) then reads 



(12) 



From Eqs. ^ and ([HI), ((A)) becomes the same as the 
thermodynamical average 



(13) 



Dynamical correlation function for an arbitrary oper- 
ator A is given as 



x(^) = — —T 



w - "H -f e„ 4- 27 



A\en), 

.... 

where L is the system size and 7 is an infinitesimally 

small energy. In our DMRG simulations, we take 7 to 
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be a small but finite value, which makes a delta function 
to a Lorentzian-type broaden function. To calculate the 
dynamical correlation function, we introduce an integral 



with 



q = 



(n=l,2,...,i) 



(23) 



1 



T^mo J- 



1 



Lj — H + e + i'-f 
(15) 

where |e) = S{£ — H)\^). The integral x'{u!) is expected to 
converge to xi'-^) after random sampling and averaging 
mentioned above.— 

During DMRG procedure before the final calculation 
of Eq. (|15p , we make use of the following approximation 
to prepare a target state for dynamical quantity: 



We note that open boundary condition is better in get- 
ting a good convergence than periodic boundary condi- 



1 



uj — 71 + e + i-f 



A\e) 



1 



uj-H + E + i^ 



(16) 

where E = (C|H|0.^ 

In dynamical DMRG calculations, we use three target 
states: |eo), ^|eo)i and [uj — H + €q + i^]^^A\eQ) at zero 
temperature, and |^), A\£^), and [uj — "H + E + ij]^^A\£^) 
at finite temperatures. 

To study the chiral order, we calculate a chirality cor- 
relation function, which is defined by 



(17) 



where S" = 1/2 and II (Ir) denotes the left (right) site of 
chirality. 



L 



II ^ — - int — 



ri 



and 



In = h int I -^-i — 

" 2 I 2 



(18) 



(19) 



with ri being the distance between both sites. The inte- 
ger part of X is represented by int(a;). 

We also examine spin-correlation function, whose mo- 
mentum representation is assumed to be 



(20) 



Dynamical chiral-correlation function is defined by 

N 



Kiq,Oj) 



n — 1 

x3(e„|K(g) 



^iiq)K), (21) 



where K(g) is the momentum representation of chirality 
under open boundary condition: 



Kiq) 



(22) 



tio: 
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IV. RESULTS 
A. Dynamical chiral-correlation function 

We calculate the dynamical chiral-correlation function 
K{q, uj) to confirm chiral properties of the ground state 
at a J = 0.5 (the chiral phase) and aj = 1.0 (the dimer 
phase). Figure [ija) shows K{q,u!) in various sizes in the 
chiral phase (a,/ — 0.5). The momentum q is taken to be 
the smallest one on each system, i.e., q = tt/{L + 1). 
The factor 7 in Eq. (PT|) . which makes a delta func- 
tion Lorentzian shape, is taken to be O.OIJ2. Since a 
peak structure in K{q,uj) is well described by a single 
Lorentzian with the half width at half maximum, 7, its 
peak position is regarded as the position of the delta func- 
tion. The peak position approaches zero energy a; = 
with increasing the system size. The size dependence of 
the peak position is plotted in Fig. [2] We find that an ex- 
cited state forming the peak degenerates with the ground 
state in the thermodynamic limit L ^ 00. This is simply 
due to the presence of degeneracy in the ground state as 
a result of chiral LRO. 

In Rcf. 10, excited states above the peak position 
have been examined by the exact diagonalization (ED) 
method and no gap has been identified between the 
ground state and the excited states. In the present cal- 
culation, it is difficult to identify precise energy of such 
excited states. For example, there is a very small hump 
structure at a; ~ 0.25J2 for L = 256 in Fig. [IJa) but 
precise determination of the energy of a corresponding 
eigenstate is impossible unlike the ED methodJ^ 

Peak positions for aj = 1.0 (the dimer phase) in 
Fig. IHb) approach a finite energy with increasing the 
system size accompanied by the reduction in spectral 
weight. Figure m shows that gap magnitude in the ther- 
modynamic limit is w ~ O.O5J2. The spectral weight at 
the gap energy seems to vanish in the thermodynamic 
limit [see Fig. [ijb)] , in contrast to the case of the chiral 
phase. 

In order to investigate the effect of temperature on 
chiral excitation, we calculate K{q, uj) at finite tempera- 
tures for the L ~ 256 zigzag chain. At aj = 0.5, spectral 
weight decreases monotonously with increasing temper- 
ature as shown in Fig. [3ja). On the contrary, K{q,Lo) 
for the dimer phase [Fig. [3]Jb)] shows a non-monotonic 
behavior: the peak position shifts to the lower energy 
side accompanied by an enhancement of spectral weight 
below Lj 0.05J2- At T/J2 = 0.04, the oj dependence is 
qualitatively similar to the case of aj = 0.5 at the same 
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FIG. 1. (Color online) Dynamical chiral correlation function 
K{q,uj) in spin-^ zigzag XY chains with various system size 
L at zero temperature, (a) aj = 0.5 and (b) qj = 1. The 
broadening factor 7 — O.OIJ2. The momentum q is the small- 
est one defined in each system. The maximum DMRG trun- 
cation number is m = 400 for L = 256. 



temperature. There is a crossover from gapped to gapless 
excitation between T/J2 =0.02 and 0.04. 



B. Static correlation functions 

In order to understand the behaviors of dynamical 
chiral correlation mentioned above, we calculate tem- 
perature dependence of static chiral-correlation function 
Cfiir) as a function of distance r. Figure [4] shows (^^(r) 
up to r = 100 for the L = 256 system. We found that 
finite-size effects seem to be serious when r is more than 
100, judging from the fact that the value of the correla- 
tion is independent of system size up to r ~ 20 and r ^ 40 
for L = 64 and L = 128, respectively. Concerning con- 
vergence of the data in terms of the truncation number 
m in LT-DMRG, we checked that m — 500 gives a good 
convergence even for the highest temperature T = 0.06 J2 
examined. At zero temperature for aj = 0.5 shown 
in Fig. IHa), Cr,(r) is almost constant, representing the 
presence of chiral LRO.— We note that the data show 
a very slight decrease beyond r ~ 80. This is actually 
the finite-size effect. In contrast to the zero-temperature 
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FIG. 2. Peak position of the spectra in Fig. [T] as a function 
of the inverse system size. 
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FIG. 3. (Color online) Temperature dependence of dynamical 
chiral correlation function K{q^Lo) in a 256-site spin-i zigzag 
XY chain for (a) aj = 0.5 and (b) a,j — 1. The broadening 
factor 7 — O.O5J2. The maximum DMRG truncation number 
is m = 500. 



result, Cnir) decreases with r at T/J2 ~ 0.02, indicat- 
ing that the LRO at zero temperature disappears at fi- 
nite temperature. The difference between the T = and 
T — O.O2J2 results becomes clear when we plot the data 
with a log- log graph (not shown) : The T = data exhibit 
an upward bending that is an indication of LRO (Rcf. 4) 
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FIG. 4. (Color online) Temperature dependence of chiral cor- 
relation function C^ir) in a 256-site spin-i zigzag XY chain 
for (a) Qj — 0.5 and (b) aj — 1.0 as a function of the distance 
r. The maximum DMRG truncation number is m = 500. 



while the T = 0.02 J2 data do not show such an upward 
bending. With further increasing temperature, Ck(^) de- 
creases monotonously. We note that the monotonous de- 
crease in C„(r) corresponds to the monotonous decrease 
in K{q, uj) with temperature as shown in Fig. [3ja). 

We also performed the LT-DMRG calculate of an order 
parameter (kx,/2) for L = 400 chain, according to a tech- 
nique proposed in Ref. 7 (not shown here). We found 
that, at T/Ji = 0.01, {kl/2) becomes zero. Together 
with the results of Ci^{r) in Fig. [DJa), we conclude that 
no chiral LRO at finite temperature in the chiral phase of 
the spin-i zigzag chain. This is consistent with the fact 
that there is no LRO at finite temperature in ID systems 
with short-range interactions. 

Figure ID^b) shows CK{r) for aj = 1.0 without chi- 
ral LRO. At T = 0, G^{r) rapidly decreases with r as 
expected. With increasing temperature, Cre(r) at all r 
increases and then decreases, showing maximum values 
at T/J2 = 0.03. This temperature corresponds to the 
crossover temperature in K{q,u}) discussed above. The 
enhancement of C^ir) can be attributed to the contribu- 
tion of excited states with strong chiral correlation. It 
is natural to assign the eigenstate at the gap position of 
K(q,uj) as such an excited state. In fact, corresponding 



FIG. 5. (Color online) Temperature dependence of spin corre- 
lation function Cs{q) in a 256-site spin-i zigzag XY chain for 
(a) a,j — 0.5 and (b) aj — 1.0 as a function of the wave num- 
ber q. The maximum DMRG truncation number is m = 500. 
Inset shows the height of the peaks as a function of tempera- 
ture. 



gap energy uj ~ 0.05 J2 is close to the crossover tempera- 
ture. Therefore, this temperature-induced chiral correla- 
tion together with the enhancement of low-energy chiral 
excitations in K{q,uj) is a demonstration of the presence 
of a chiral state in excited states in the dimer phase. 

It is known that chiral LRO at the ground state does 
not need LRO of spin^ In fact, spin-correlation function 
decays in the power law while chiral-correlation func- 
tion remains finite in the limit of r — >■ 00. However, 
the magnitude of the spin correlation seems to be cor- 
related with that of the chiral correlation)^ Therefore, it 
is expected that the spin correlation shows temperature 
dependence similar to that of the chiral correlation. In 
order to confirm this, we also calculate spin correlation 
function Cs{q) in Eq. ((20)) . 

Figure [5] shows Cs{q) for aj = 0.5 and 1.0. For 
aj = 0.5, peak height at g ~ 0.537r decreases with in- 
creasing temperature. The decrease is similar to that in 
the chiral correlation shown in Fig. HJa). In the case of 
the dimer phase {aj = 1), peak height increases up to 
the characteristic temperature T/J2 ^ 0.03 before it de- 
creases, similar to the behavior of C^ir)- The correspon- 
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dence between the chiral and spin correlations, thus, con- 
firms the presence of excited states that have strong chiral 
correlation in the dimer phase. The position of the peaks 
is not the same for aj = 0.5 and for aj = 1. The posi- 
tion is consistent with a characteristic vector of the spiral 
phase in the classical limit, given hy q = cos^^(— aj/4). 
We also note that the peak at q = and dip at g = tt 
for a J = 1.0 seen in Fig. [Hb) can be assigned to the 
signature of dimer correlation. 

The most important finding in this section is the fact 
that the chiral correlations arise at finite temperature 
from the dimer ground-state phase. This reminds us 
of similar behaviors seen in order-by-disorder phenom- 
ena, i.e., induced ordering at finite temperature from dis- 
ordered ground state.— However, physics behind the 
present results and order by disorder looks different, be- 
cause no disorder exists in the ground state in the present 
case. 



V. CONCLUSION 

We have investigated finite-temperature effects on vec- 
tor chirality in spin-i zigzag XY chains, employing LT- 
DMRG. We have compared two cases: with and without 
chiral LRO. The latter corresponds to the dimer phase. 

We have calculated the temperature dependence of dy- 
namical and static correlations of vector chirality as well 
as static spin correlation. In the chiral phase, chiral 
LRO disappears at finite temperature as expected form 



the fact that there is no LRO at finite temperature in 
ID systems with short-range interactions. In the dimer 
phase next to the chiral phase, static chiral correlation 
as well as spin correlation is enhanced with increasing 
temperature up to a characteristic temperature that is 
related to the gap magnitude of chiral excitation. This 
temperature-induced chiral correlation is a demonstra- 
tion of the presence of a chiral state in excited states. 

The present result also reveals usefulness of LT-DMRG 
for the numerical studies of ID frustrated quantum spin 
systems at low temperature. 
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